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Using the Dyson equation of motion for both the diagonal one-particle Green function, G(k, lS) 
and off-diagonal Green function, F(k,tu), at the level of the Hubbard-I decoupling scheme, we have 
found that they have four poles symmetric in pairs, justifying a more elaborated calculation done 

■ by the Zurich group by means of the T-Matrix approach (Pedersen et al, Z. Physik B 103, 21 
(1997)) and the moment approach of Nolting (Z. Physik 255, 25 (1972)). We find that the energy 
spectra and the weights of G(k,u) and F(k,u>) have to be calculated self-consistently. G(k,u>) sat- 

D ' isfies the first two moments while F(k, ui) the first sum rule. Our order parameter a(T) is given by 

5^ , l/Ns ^2^e(k)A(k). Due to the fact that we have a purely local attractive interaction A(fe) can be 

■ of any s-type wave. However, for a pure s-wave, for which a(T) = 0, we go back to the mean-field 
BCS results, with a renormalized chemical potential. In this case, the off-diagonal Green function, 
F(k, ui), satisfies the first two off-diagonal sum rules. We explicitly state the range of validity of our 
approximation . 

. Pacs numbers: 74.20.-Fg, 74.10.-z, 74.60.-w, 74.72.-h 

o 

o . 

After the discovery of the high-T c materials the study of correlations has gained interested due to the fact 
| that there is the belief Q that the normal properties of these materials could be explained in the framework of the 
Hubbard model |,|, since electron correlations are strong, i.e., the on-site electron-electron repulsions U are much 
larger than the energies associated with the hybridization of atomic orbitals belonging to different atoms j5j. This 
Hamiltonian is a kind of minimum model || which takes into account quantum mechanical motion of electrons in a 
solid, and nonlinear repulsion between electrons. Even though this model is too simple to describe solids faithfully, 
serious theoretical studies have revealed that to understand the various properties of it is a very difficult task. Its 
qq ■ study will prove useful in developing various notions and techniques in statistical physics of many particle physics. 
^\ \ Since the high-T c superconducting materials are extreme type II superconductors, with a short coherence volume, 

■ one might take this fact as an indication of tightly bound pairs or/and that correlations effects strongly affect the 
properties of such materials. However, this scenario has been challanged by recent tunneling measurements |t| which 
lead to the conclusion that, for example, underdoped i^SV^CaC^Os-iS is described by intermediate coupling 
interaction, because the pairing fluctuations persist up to T* where T* (T* > T c ) is the temperature of pair 
fluctuations and T c is the superconducting critical temperature. 

One of the simplest model featuring superconductivity and allowing a systematic study of correlations is the 
attractive Hubbard model || which we adopt in this paper as the counterpart of the usual Hubbard model. This 
model has been used to explore qualitative features of the superconducting phase transition . Ref . p3| is mainly a 
review of the analytical work done on this model. Recently, Huscroft and Scalettar [[yj find that the superconducting 
order parameter is more stable that the charge density wave order at half-filling in the presence of disorder. Then, 
due to these considerations, we concentrate in the superconducting properties leaving outside any treatment of charge 
^ , density wave order. 

We will use the Dyson equation of motion technique for both the diagonal, G{k,ui), and off-diagonal one- 
particle Green function, F(k,u>), at the level of Hubbard-I decoupling scheme j|. The main theoretical conclusion 
of this paper is that both the diagonal and off diagonal one-particle Green functions have four poles. These poles are 
symmetric in pairs verifying a more elaborated calculation of the Zurich group |l3| . 

The model we study is the Hubbard model || 

H = tr-rd c-> H n?- — acL o , (1) 
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where cl (c^) are creation (annihilation) electron operators with spin a. n- i(y = ct c* a . U = —\U\ is the local 
attractive interaction and p, the chemical potential (we work in the grand canonical ensemble). We have adopted 
Einstein convention for repeated indices, i.e., for the N s sites i, the z nearest-neighbor sites j and for spin up and 
down (<t = —a = ±1). t^j = —t, for n.n. and zero otherwise. Other types of hopping, i.e., t' ^ between next nearest 
neighboors {n.n.n.) could be considered in our formalism |l4j ]. In this paper, we restrict ourselves to n.n. hoping. 
We need to evaluate the equation of motion Q for the operators cj CT and ct . They are 

dc- del r + 1 

where the sign — means the commutator. Combining Eqs. (0,0) we 6 e t 



Next, the one-particle Green's function is defined as 

G ff (?, J ; r) = - ( (T T c ?ct (t) ; c* (0) ) ) , (4) 



where XV means time ordering. Combining Eqs. (j^H), and Fourier analyzing the time and space variables we end up 
with the following equation for G(k,u>) 

(uj-e j :)G(k,uj) = l + Ur {2) (k,oj) , (5) 

where T^(k,Lu) is the Fourier transform of the doubly occupied Green function [ JDS] 

r( 2 )(fc,c) S ((r^(r) C?CT (r); C t CT (0)))(fc, W ) . (6) 

As we see from Eq. (^), the one-particle Green function and the doubly occupied Green function are connected 
thruout the equation of motion. We get the doubly occupied Green function mainly due to the presence of four 
operators in the Hubbard interaction. In Eq. (^), n is = c\ s c is is the occupation number operator. Our next step is 
to apply the Dyson equation to T^{k,uj), which is given by 

(w - U)T^{k^) = P» + * ?) r((n & (T)c r<r (r);^ <r (0))) (Z,u>) , (7) 



which we have obtained under the assumption that 

On? 



l. a 



dr 



. (8) 



Eq. (||) is certainly an approximation which reproduces the Hubbard-I solution in the equation of motion approach 
and it is only valid at the level of Eq. (Q) . We leave for the future |l(J the study of the effect of going beyond the 
approximation given by Eq. (||) . Now we perform a decoupling in the spirit of Hubbard Q] as follows 

((^(r)c r >); c+ ct (0))) » p 9 G Q (r) - . (9) 

Combining Eqs. ([|||), we arrive to the following expression for T 2 (k, u>) 

( LU -U)T^(k,Lo)=p rf + p^G(k,Lu)-^lF^k,cu) , (10) 



where a{T) = l/N s J2]:s(k)A(k). Combining Eqs. (|,QQ we get 

[(u>-E i )(u-U)-p 9 Ue i ]G(k,w)=w-U(l-p 9 )-afflF*(k,w) . (11) 

Thus, from Eqs. (|TT|) , we see that G and F are coupled and that G reduces to the Hubbard-I solution, as it should 
be, when a(T) = 0. Let us pause for a while to explain the notation. The parameter a(T) is going to be our order 
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parameter, in complete analogy with the decoupling scheme in mean-field treatments (BCS one, for example). Next, 
we have to find the time evolution for F^(k,w). A similar analysis, i.e., another Hubbard-I decoupling scheme for 
F'(k, cj), shows that 

[{u + eg) {u + U)~ PaUej:] [k, u) = a* (T)G(k, to) (12) 
Eqs. (0,[l2|) produce for G(k, u>) and F^(k, w) the following solutions, 

[lu - U{1 - pa)) [(lj + e % ){u + U)- P aUe % ] 



G(k,uj) = 



[(w - e s )(u - U) - pallet [(w + e £ )(u + U)- p s Ue t ] + \ a(T) | 2 



Fi k co) = a{T) [U U{1 Ps)] (13) 

V ' ; [{w-e^-U)-p 9 Ue^[{u} + e i ){y} + U)-p 9 Ue^ + \a{T)\^ ' V ' 

From Eqs. (|l|) we conclude that, for a(T) 7^ 0, 

G (*»=E^#k • F *">=S^fk • (14) 

i.e., the one-particle Green functions have four poles (lifetime effects are neglected here) which turn out to be 

&i(k) = - Cl 2 (k) = uj (k) , Cl 3 (k) = - Q 4 (k) = ui(k) , (15) 



where 

<i(fc) = \ 



1/2' 



C(k) ± [c 2 (k) - 4 (| a(T) | 2 +(1 - Ps ?U 2 ef) 
C(k) = U l + e\ + 2Ue l ; e % = e{k)-p , (16) 

with e(k) = — 2tY2j = i cos(kj) and d the lattice dimension. Taking a closer look to w^(fc) in Eqs. ( |l6| ) we conclude 

that these poles have almost the form of the poles for G(k,u) and F(k,u>) obtained in Ref. |l3|, since they give 
four solutions, symmetric in pairs, respecting what we call the BCS symmetry. In Ref. |l3f| , we interpreted these 
two symmetric solutions as corresponding to the BCS solution (the opening of the BCS gap around the chemical 
potential) and to the pair physics (the correlation gap), respectively. In other words, our two Hubbard-I decouplings 
have given an additional contribution, which is due to the presence of pair fluctuations above T c and which remain 
for T < T c . These pair fluctuations are the electrons which are not in the Meissner state. The only qualitative 
difference with respect to the results of Ref. ]T3] is that here we do not have lifetime effects. 
The spectral weights (Eqs. (|l4|)) are given by 

&i(k) 
a 2 {k) 



[uJo(k) 
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and 



2 Wl (fc)( W 2 (fc)-^(fc)) 



h(k) = ^T)Mk)- U(l-p 9 )] _ a{T) - {%) 
2co (k)(^ (k) - L0\{k)) 



P 2 (k) = ^T)Mk) + U(l-p 9 )] _ a{T) - {%) 
2u (k)(cj (k) - wf(fe)) 
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2w x {k){<jl{k)-u>l{k)) 

2iOx{k)(ul{k) - ul{k)) 
Now it is an easy matter to convince ourselves that the following relations are satisfied 

di(ft) + d 2 (ft) + d 3 (fc) + d 4 (fc) = 1 , 
w (ft)(di(ft) - da(ft)) + ui(k)(a 3 (k) - 64 (ft)) = + p s U , 
0i(k) + fo(k) + p 3 (k) + fa(k) = , 



(18) 



(19) 



which are the first three sum rules for the moments ]T^,0 . The second off-diagonal moment is not satisfied because 
our order parameter is not A(T) but a(T). A discussion of the failure of a Hubbard I-type solution (for a(T) = 0), 
as the one presented in Eqs. (|ll|), has been pointed out by Laura Roth Q many years ago, where she remarked that 
the diagonal Hubbard I solution only satisfied the first two sum rules for the moments. We should point out that the 
moment solution and the Hubbard I-type solution are different approaches, which turn out to be approximately equal 
due to the fact that some of the first sum rules are satisfied. Laura Roth's criticism to the Hubbard-I solution is still 
valid in the present calculation since we should have a correlation gap for any finite value of U ft. In order to close 
this gap, lifetime effects are called for. 
From the spectral theorems we have |l9| 



-f-oc 



Aij(uj)du> 
exp(/3u)) + 1 



( c n c 3i> 



Bi j (uj)duj 
exp(/3uj) + 1 



Then, we have 



-E 



A(k, lu)cLj 
exp(/3oj) + 1 



a(T) 



-E 



-Hoc 



e(k)B(k,uj)duj 
exp(/3uj) + 1 



(20) 



(21) 



where (5 — l/T is the inverse of the temperature and a(T) has been defined just after Eq. (jl(]). A(k,uj) and B(k,iv) 
are the one-particle spectral densities given as 



A(k,cj) = — lim Im[G(\s.,oj + iS)] ; B(k,oj) = lira Im[F(k, oj + iS)] 

7T 8^0+ ' 7T 8^0+ 



(22) 



In consequence, by combining Eqs. ( pO| , pT| , p2[ ) with our Green functions (Eqs. (|l3|)), we obtain the following 
self-consistent equations 



1 

U 



+4 



eJV(e) 



-4 
-4 



N(e) 



/3i(e) + /3 2 {e)ex P (/3u {e)) f3 3 (e) + P^exp^uj^e)) 

exp((3u) (e)) + 1 ' exp(/3u>i(e)) + 1 
&i(e) -I- a.2{e)exp{(3uj {e)) 613(e) + a i {e)exp((3uji(e)) 



exp((3uj (£)) + 1 



exp((3uji(e)) + 1 



ds 
de 



(23) 



where N(e) is the 2D density of states, t = 1, and we have chosen p = p a = p 9 , i.e., we are in the paramagnetic phase. 

We indicate that our Eqs. (23) will respect particle- hole symmetry ^(J. From the analysis of the first equation 
of Eq. d25| ) and the definition of the our order parameter we can conclude that it allows any type of s-type of wave 
symmetry, since when performing the ^-integration of A(fc) x e(k), we see that a(T) ^ only if A(fc) is of s-type. 
So, in this case, Eq. (^3|) can give rise to an order parameter of symmetry different from pure s-wave, a conclusion 
which was reached in a previous work plj using the sum rules both for the diagonal and the off-diagonal one-particle 
spectral functions. However, for a pure s-wave, i.e., A(fe) = const., we get a(T) = 0. So, our approximation fails psj] 
and we must go back to Eq. (Js[) and directly perform our approximation in T%(k, u>). We arrive to the mean-field 
BCS results, where the chemical potential gets renormalized by the Hartree shift, i.e., pU . Due to our lazyness of 
languaje, we have used the word order parameter in this paragraph to denote A(fc), even though we defined a(T) just 
after Eq. ( |Tl] ) as the order parameter of our theory. Naturally, a(T) is the integral of A (ft) weighted with e(fc). So, 
they are related in some way, i.e., if A(fc) = 0, then a(T) = 0. 
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In short, using the Hubbard-I decoupling scheme for both the diagonal and off-diagonal one-particle Green functions 
we have shown that these Green functions have four poles, symmetric in pairs, which qualitatively verify the more 
elaborated calculation of Ref. Q|, as it has been previously discussed. Our one-particle Green functions satisfy sum 
rules for the moments and we have obtained other symmetries than a pure s-wave order parameter. The range of 
validity of our approximation is contained in Eqs. ||). Now we solve our Eqs. ( p3| ) in a low order approximation: 



We fix the value of our order parameter, a{T), and find the chemical potential using the first of Eqs. (23). For 



p = 0.01, U/t = -8.0 and ^r 1 = 0.1, we nnd A«A ~ -2.825. We should say that our approach is valid for \U\ > W, 
where W — 8t is the bandwidth in two dimensions. Our calculation have been performed for /3 = 1/T = 100.0. We 
leave for the future ]l6| the numerical evaluation of the critical temperature, T c , and the order parameter as function 
of temperature for different values of U/t and electron concentration. 

In order to deal with ci-wave superconductivity, we should study a nearest neighboor (n.n.) attractive interaction. 
This model has been considered few years ago, at the mean field level, by Meintrup, Schneider and Beck |^3) (See, 
also Ref. We mention, while leaving, that lifetime effects can be included in a natural way both in G(k,ui) and 

F(k,uf), as it has been previously done in Ref. ]2j| for the self-energy. Work along these lines is in progress. 
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